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Hard-sphere colloids are popular as models for testing fundamental theories in condensed matter 
and statistical physics, from crystal nucleation to the glass transition. A single parameter, the 
volume fraction (4>), characterizes an ideal, monodisperse hard-sphere suspension. In comparing 
experiments with theories and simulation, researchers to date have paid little attention to likely 
uncertainties in experimentally-quoted (f> values. We critically review the experimental measurement 
of <j> in hard-sphere colloids, and show that while statistical uncertainties in comparing relative values 
of cf> can be as low as 10~ 4 , systematic errors of 3-6% are probably unavoidable. The consequences 
of this are illustrated by way of a case study comparing literature data sets on hard-sphere viscosity 
and diffusion. 
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I. INTRODUCTION 



Excluded volume effects dominate the behaviour of 
liquids around the triple point [l[ and play a key role 
in structuring crystalline [2| and amorphous solids Q. 
Thus, hard spheres have long functioned as a reference 
system for theoretical and simulational studies of con- 
densed matter. In 1986, Pusey and van Megen demon- 
strated HI that suspensions of sterically-stabilised poly- 
methylmethacrylate (PMMA) colloids showed nearly- 
perfect hard-sphere equilibrium phase behaviour, under- 
going a first-order phase transition from a fluid to a crys- 
talline state at concentrations around those predicted 
some time ago by computer simulations for hard spher- 
ical particles [H, Q . Soon afterward, the same authors 
showed Q that PMMA colloids underwent a glass tran- 
sition at even higher concentrations. Subsequently, in 
a 1991 review, Pusey enunciated the 'colloids as atoms' 
paradigm — Brownian suspensions can be used as 
'test tube simulations' of many generic condensed mat- 
ter phenomena such as crystallization and vitrification. 
Since then, the use of colloids as models has become very 
popular, especially since the addition of non-adsorbing 
polymers can be used to induce an inter-particle attrac- 
tion 'tuneable' separately in its range and depth [{J [lj| ■ 
In such colloid-polymer mixtures, gas-liquid coexistence 
can be studied [III EH > including interfaces and critical- 
ly [HI 13: as wen as novel modes of arrest [Hj]. Due 
to their slow intrinsic time scales, colloids are also ideal 
models for studying phase transition kinetics (l6| . 

While a range of inter-particle interactions are now 
available in model colloids, hard spheres remain an im- 
portant reference system for which very direct compari- 
son between experiments and theoretical calculations or 
computer simulations is in principle possible. The be- 
haviour of a single-sized, or monodisperse, system of hard 



spheres is controlled by one parameter, the volume frac- 
tion (j>, i.e. the fraction of the total volume V that is 
filled by N spheres, each of radius a, 
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Since <f> is precisely known in theory or simulations, a 
comparison with experiments is straightforward provided 
that this quantity is also accurately measurable for real 
suspensions. Much of the literature has indeed pro- 
ceeded on this basis, assuming that <f> is unproblemati- 
cally known from experiments. 

However, as Pusey and van Megen pointed out in a 
symposium article 17 following their Nature paper, Q 
the experimental determination of <f> is emphatically not 
unproblematic, because: (1) no real colloid is truly 
'hard', since there is always some softness in the inter- 
particle potential; and (2) real colloids always have a 
finite size distribution, i.e. they are polydisperse. Thus, 
Pusey and van Megen calculated an experimental 'ef- 
fective' hard-sphere volume fraction </>e, and found that 
the freezing and melting volume fractions of their system 
were 0.494 and 0.535 respectively, compared to 0.494 and 
0.545 in simulations.]^ [f| Their careful conclusion reads: 
'Despite ambiguities ... in the experimental determina- 
tion of the coexistence region this difference is probably 
significant.' 

The same degree of caution has not characterized the 
literature since. Experimental reports typically do not 
discuss in any detail the method used for arriving at <f>. 
On the other hand, theory or simulations almost always 
take experimental reports of <fi at face value and proceed 
to use the data on this basis. This is unsatisfactory, par- 
ticularly in situations where theory testing demands a 
degree of accuracy and certainty in the experimental <j> 
that is probably unattainable. In this paper, we critically 
review a plethora of methods for the experimental deter- 



mination of 4> m hard-sphere suspensions, evaluate the 
degree of accuracy attainable in each case, comment on 
the potential discrepancies between methods, and give 
study showing how different experiments should 
be compared taking into account possible difference in <f> 
determination. 

With the increasing popularity of confocal microscopy, 
direct counting of particles is becoming a standard 
method for determining </> (see Sec. IV D[) . This method 
depends on knowing the particle size. Thus, after in- 
troducing model colloids (Section HI]) . we review particle 
sizing (Section The 'classic' method for determin- 

ing (j) is via the crystallization phase behaviour, which 
changes with polvdispersitv[l8l fl9| . So we review poly- 
dispersity measurements (Section IIV[) before turning to 
consider the determination of <j> in detail (Section|V|. We 
finish with a case study (Section IVI[) and a Conclusion. 



II. THE PARTICLES 

We focus on suspensions of nearly-perfect hard spheres. 
A system that potentially behaves most like perfect 
hard spheres is charge-stabilised colloidal silica in wa- 
ter. When the char ges are sufficiently screened out by 
the addition of salt [20l - |22| . the resulting suspension is 
very close to being hard-sphcrc-like. However, a signifi- 
cant drawback of silica as a model system is that these 
particles have a density p rts 2.2 g/cm 3 , and it has proved 
impossible to find solvents that match this density. The 
sedimentation problem can be alleviated by using smaller 
particles, but at the expense of increasing polydisper- 
sity. Charged-stabilised polystyrene spheres are also po- 
tentially very hard-sphcre-like, and are easier to density 
match, though (unlike silica) almost impossible to refrac- 
tive index match. 

A more popular model hard sphere system is PMMA, 
sterically stabilized by a 6 5, 10 nm layer of PHSA (poly- 
12-hydroxystearic acid) [2j|. This layer confers a degree 
of softness to the interparticle potential on the scale of 
S. PMMA particles can be dispersed in solvent mixtures 
that match both the particles' density and index of re- 
fraction [24], [2f| . Unfortunately, particle swelling by sol- 
vents is endemic [H, [2?| ■ Furthermore, the swelling pro- 
cess can take several weeks, so the particle size changes 
over time, though heat shock may speed up the process 
to taking only a few hours [28j . As swelling is poorly 
characterized, in situ measurement is the only reliable 
means of ensuring that it is complete before the parti- 
cles are used in experiments. This is particularly impor- 
tant at high (j), where many properties are steep func- 
tions of the concentration: an x% increase in the particle 
radius translates into > 2>x% in <f). Thus, e.g., an index- 
matching mixture of cis-decalin and tetrachloroethylene 
causes 20% swelling, which has a drastic effect upon the 
colloid volume fraction [27} ■ Finally, batch-to-batch vari- 
ations generate further uncertainties. 

Below, we focus on PMMA particles, although much 
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FIG. 1: Definitions of particle radii, (a) Non-sterically- 
stabilized particles, e.g. silica, with a core radius a c . (b) 
Sterically-stabilized particle, with surface 'hairs' (not drawn 
to scale), where additionally the average hair thickness 5 and 
the core-shall radius acs ~ ac + S axe needed for a full char- 
acterisation. 

of what we say will also apply to silica and other model 
systems. 

III. MEASURING SIZE 

The basic parameter characterizing spherical colloidal 
particles is their radius, a. Some methods of determining 
</>, most obviously by 'counting' from confocal microscopy 
images, depend directly on measuring a. In this section, 
we critically review the measurement of particle size. 

A. Definition of particle radius 

At first sight, defining 'the particle radius' should be 
simple in a nearly-monodisperse colloid. But this is de- 
ceptive, Fig. Q] For sterically-stabilized particles like 
PMMA, we can usefully define as least four different 
radii. First, the hydrodynamic radius, an, occurs in the 
relationship between the drag, /, on a particle moving 
at velocity v in a fluid of viscosity rj at low Reynolds 
numbers, f = £v, where the friction coefficient 

£ = 6irr]a H . (2) 

Relating an to is a non-trivial problem in hydrody- 
namics. Somewhat more directly related to <f> is the core 
radius, ac, which is the radius of the sterically-stabilized 
particles minus the stabilizing hairs. Thirdly, if we can 
determine the average hair thickness 5, then the core- 
shell radius, acs = ac + S. For 'hairless' particles, such 
as charge-stabilized silica, acs = o-c- Finally, we may 
assign an effective hard sphere radius, a e ffi to the par- 
ticles to obtain the best fit to theory or simulations of 
hard-sphere behaviour in a certain range of cf>, so that 
a e ff is inevitably dependent on the chosen property and 
<p range. 

We will not review established sizing methods in any 
detail, but will reference existing literature and note cau- 
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tionary points. Then wc will introduce a number of newer 
methods. 



B. Measuring radius: established methods 

Scatterin g m ethods have a long history in sizing spheri- 
cal particles }29l - [3l| . Static and dynamic scattering deter- 
mine the size of particles by measuring the time-averaged 
or fluctuating intensity of the scattered light respectively. 

Dynamic light scattering (DLS) and its X ray equiv- 
alent, X ray photon correlation spectroscopy (XPCS), 
measure the diffusion coefficient of particles, which is re- 
lated to the friction coefficient via the Stokes-Einstein- 
SutherlandrelationHHH: D = k B T/£. DLS and XPCS 
therefore determine an (cf. Eqn. [2]), and are most useful 
in the case of particles consisting of core only, such as sil- 
ica, since it is less clear how to relate an to acs for core- 
shell particles such as PMMA. Note that the accuracy 
of this method depends on having an accurate value for 
rj, the solvent viscosity, which is temperature dependent. 
For example, we have found that for the common sol- 
vent mixture cyclohexylbromidc and decalin (85%/15% 
by weight) 24°C, 77 = 2.120 mPa-s and d-q/dT = -0.029 
mPa-s/K. Thus, a 1°C uncertainly in T is a 0.3% uncer- 
tainty of T but a 1.7% uncertainty in T/rj and therefore 
in an- 

Static light scattering (SLS), small-angle X ray scat- 
tering (SAXS) or small-angle neutron scattering (SANS) 
can potentially determine ac and acs- Since the core 
and shell of (say) a PMMA particle in general has differ- 
ent contrasts to light, X rays (refractive index, n, in both 
cases) and neutrons (scattering length, b) , the diffraction 
pattern of a single particle is determined by the interfer- 
ence of radiation scattered from these two parts. Fitting 
this diffraction pattern (the form factor) therefore can 
in principle yield ac and acs = a c + S. In a solvent 
with n or b quite different from both the core and the 
shell, the whole entity scatters more or less as a homoge- 
neous sphere and a radius close to acs is returned from 
form-factor fitting. When solvent mixtures are used to 
'tune' the relative contrasts of core and shell, even a small 
amount of a minority component in the solvent mixture 
can swell the particles by up to 10% or more [26l [271] . and 
the fractional swelling of core and shell is not necessarily 
identical. In XPCS, where the shell has little contrast, 
S cannot be accurately determined; however, the bright- 
ness of the beam gives many orders of oscillations in the 
form factor, allowing very accurate data fitting. 

For both static and dynamic scattering, samples must 
be dilute enough so that the properties of non-interacting 
particles are measured in the single-scattering limit. The 
only sure way to know that this has been achieved is to 
collect data at different and look for the convergence 
in the <j) —> limit. Static scattering at finite </> gives the 
static structure factor as a function of scattering vector, 
S(q). Fitting this to, e.g., the Percus-Yevick form (3f| 
or simulations yields simultaneously a c g and 0, although 



polydispersity is a significant complication[36j . Alterna- 
tively, the Bragg peaks in S(q) from colloidal crystals at 
fluid-crystal coexistence can be used to deduce a e s if the 
melting point is known (but see Section ITVl for caveats). 

Electron microscopy (EM) measures ac of dried parti- 
cles, because drying collapses the steric-stabilizing 'hairs' 
in core-shell particles such as PMMA, and deswells par- 
ticles swollen by solvent when dispersed. Various optical 
microscopies can, in principle, be used in the same way 
as EM for sizing particles; caveats are pointed out in 
Section HIDJ 

C. Measuring radius: newer methods 

1. Differential dynamic microscopy 

DLS measures diffusion via determining the intermedi- 
ate scattering function (ISF) , which is the spatial Fourier 
transform of a time-dependent density-density correla- 
tion function [29l - [3l| ; it requires a laser and bespoke elec- 
tronics (a correlator). Recently, a method for measuring 
the ISF has been demonstrated [37| that requires only the 
use of everyday laboratory equipment, viz., a white-light 
optical microscope and a CCD camera. This method, dif- 
ferential dynamic microscopy (DDM), exploits the fact 
that the intensity of a low-resolution microscope imag- 
ing is linearly related to the density of particles in the 
sample being imaged [38j. Thus, correlating the Fourier 
transform of the images gives directly the ISF. 

2. Particle tracking 

Being a scattering method, DLS works in reciprocal 
space. DDM uses microscope images, but also yields the 
ISF in reciprocal space. In both cases, the measured 
quantity is the diffusion coefficient, which controls the 
mean-squared displacement (MSD) of Brownian parti- 
cles: (Ar 2 (r)) = uDt, with n = 2,4 or 6 in 1, 2 or 3 
dimensions respectively. Direct real-space methods for 
measuring the MSD are increasingly popular. The mo- 
tion of particles in a dilute sample (0 <C 0.01) can be cap- 
tured by video microscopy [Hj], and the particle motion 
tracked [40j using publicly available software (4l). Pro- 
vided the microscope has been properly calibrated, such 
tracks yield the MSD. 

Problems can occur at short and long times. The issue 
at short times is measurement error due to pixellation. 
The pixellation error for a particle whose image is N 
pixels in diameter with individual pixels of width M is 
roughly M/N. In 1 dimension, a positional uncertainty 
of X generates an apparent MSD of X 2 , which is time 
independent, giving 

(rL aB (T))=(ri uc (T))+X\ (3) 

If the short-time MSD plateau due to X 2 is observed, it 
provides an excellent means of determining the measure- 
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FIG. 2: Mean square displacement of ac = 1.5 /im 
polystyrene spheres in a water-glycerol mixture, (o): raw 
data; (+): {Ax 2 )—X 2 where the estimated noise level (dashed 
line) is -X" = 0.08 /im. Solid line: a linear fit to the raw data 
in the range 10° < r < 10 1 s, giving D = 0.0254 ^m 2 /s. 
With the noise subtracted off, a linear fit to all of the data 
at r < 10 1 s gives D = 0.0248 (im 2 /s. In this experiment, 
N = 9 and M = 0.64 /im/pixel, so the estimated noise level 
M/N = 0.07 fim is comparable to the observed X = 0.08 /im. 



determination of g(r). Finally, note that the scanning 
mechanism drifts, so that calibration on the day of mea- 
surement is important. 



4- Holographic microscopy 

A collimated laser beam directed through a microscope 
objective scatters off a particle. The scattered and un- 
scattered beams interfere in the focal plane to form a 
hologram. At low enough <fi, the holographic image of a 
single, optically homogeneous particle can be fitted using 
Lorenz-Mie theor y to determine its position, size, and re- 
fractive index [4a-l48|. This has been demonstrated with 
~ 100 nm to 10 /im particles. The radius ac of an indi- 
vidual particle can be measured to ±10 — 30 nm from a 
single snapsho t 14911 . Multiple measurements further im- 
prove on this j48[. The sizing of core-shell particles has 
not yet been attempted. Note that while this method 
fails for exactly index-matched particles, a mismatch of 
as little as 1% is sufficient to render it usable (D. Grier, 
personal communication) . 



ment uncertainty X. Figure [5] shows that it is important 
to take this term into account for accurate determination 
of D by tracking. Note that changing the parameters 
used to identify particle positions can often influence X , 
for better or worse [4p| ] . 

At long times, the measured MSD can become non- 
linear due to particles disappearing from the field of view, 
either because they leave laterally or because they be- 
come defocussed. Since the MSD at any r is computed 
based only on particles which have been observed for at 
least as long as t, too few particles may contribute at 
large r for proper averaging. 

3. Confocal microscopy 

By using a pinhole to reject out-of-focus light, laser 
confocal microscopy is capable of generating images deep 
inside (~ 100/xm) a concentrated suspension of fluores- 
cent colloids. Thus, by processing images taken scan- 
ning through a sample, a 3 dimensional image of many 
thousands of particles can be reconstructed and their co- 
ordinates obtained fill |42| . In a sample where particles 
are touching, the peak of the calculated radial distribu- 
tion function, <?(r), gives acs- Touc hing particles can 
be generated in a spun-down sediment |43|. or by induc- 
ing a very short range attraction (e.g. by adding small 
non-adsorbing polymers). While the first peak of g(r) 
is typically averaged over > 10 8 correlations, and so is 
in principle highly accurate, the particle size in a typi- 
cal confocal image is < 10 pixels, which limits the best 
accuracy to ~ 0.1 pixel. Furthermore, microscopy is sub- 
ject to certain systematic errors 0, Hj| that affect the 



D. Measuring radius: case study 

To illustrate the difficulties in pinning down a value for 
the radius of a batch of particles, we reproduce data for 15 
preparations of fluorescent PMMA colloids by Bosma et 
al.[26|. 'Wet' particles (suspended in hexane) were char- 
acterized by SLS (far from index matching, so that acs 
is measured), while dried particles were sized by EM. 
In some cases, direct measurement of sizes from confo- 
cal micrographs was performed. The peak of g(r) from 
confocal microscopy of a close-packed sample was also 
reported for one sample. We take the EM radius, always 
the smallest, to be ac- In Fig. [31 we plot A = ax — ac, 
where ax is the radius from method X. The likely length 
of PHSA 'hairs' is 6±1 nm [HoL l5li . although oligomers of 
up to 15-20 nm may be present [5l|. The hairs therefore 
account for the lower bound of A > 10 nm. Larger A 
values are likely due to particles swelling in hexane, with 
larger particles swelling more (so that in most cases the 
swelling is ~ 4% of acs)- 

The direct confocal measurements are consistently 
higher than the SLS data by > 10 nm. This illustrates 
the difficulty of direct measurements from any optical im- 
age: the image of a single particle is far from sharp at the 
edges, both due to geometric and diffraction effects. The 
measurement from g(r) is likely more accurate, since it 
relies on locating particle centres rather than edges. It 
is not clear why the one example of such measurements 
shown in Fig. [3] is also significantly higher than the SLS 
result. Overall, these data illustrates that particles sizes 
quoted in the experimental literature may be subjected 
to significant systematic uncertainties that are often not 
always reflected in the (statistical) error bars. This must 
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FIG. 3: Radii of different batches of fluorescent PMMA par- 
ticles determined by various methods [2rj|: • static light scat- 
tering, A confocal microscopy, ■ g(r) peak. 



be taken into account if the particle size is then used in 
calculating (f>. 



IV. MEASURING POLYDISPERSITY 

Polydispersity in general refers to the existence of a 
distribution of particle properties, such as size, shape, 
charge, magnetic moment, etc. Hard sphere colloids have 
a distribution of radii, P(a), for which we define the poly- 
dispersity, er, as the standard deviation of this distribu- 
tion divided by the mean: 



(4) 



Very monodisperse PMMA has polydispersities ap- 
proaching 3%, however 5-6% is typical [26j for 'monodis- 
perse' PMMA. Note that some particles, including 
PMMA, frequently display a bimodal distribution due to 
secondary nuclcation, so that a full distribution is needed 
to characterize them. 

Polydispersity is relevant here because it affects the 
equilibrium phase diagram. Monodisperse hard spheres 
freeze at <f>F = 0.494 to form crystals at the melt- 
ing point 4>m = 0.545. These two values are often 
used as fixed points for determining <p in experiments 
(see Section IV Bl) . Theory [ljl, [H] and simula- 
tions [l9l l52l l53l] show that even small a may shift (f>M 
and/or <j)p significantly, Fig. |4j Indeed, particles with 
<7 higher than some terminal value a* will fail to crys- 
tallize at all experimentally or in simulations, although 
theory fl8j predicts phase separation into coexisting solid 
phases. Simulations |54f pre dict that a* « 7%, consistent 
with early experiments [55J. Determining polydispersity 
is therefore important for measuring 4>. Note that the 
whole distribution and not just its variance may matter, 
e.g. in determination nucleation rates [56l|. 

All the methods reviewed in the last section can po- 
tentially yield information on the size distribution. Di- 
rect imaging, EM or optical microscopy, can estimate the 



FIG. 4: The theoretical phase diagram of hard spheres at 
different polydispersities, a. F = fluid, S = (crystalline) 
solid; thus FSS denotes fluid-solid-solid coexistence. Replot- 
ted fromdl. 



full P(a), subject to the same caveats already discussed. 
Moreover, larger particles may swell more (Fig. [3J see 
also [lU ) , giving a correlation between size and shrinkage 
upon drying, so that wet and dry P(a) may be different. 

In DLS or XPCS, the ISF from a hypothetical 
monodisperse suspension decays exponentially with time. 
Polydispersity turns the ISF into a sum of exponentials. 
In static scattering, monodisperse particles give sharp 
minima in the form factor, which are smeared out by 
polydispersity. (Note that multiple scattering has the 
same effect, and so can masquerade as polydispersity.) In 
principle, these features can be fitted to yield P(a)[29|. 
subject to all the usual problems and uncertainties as- 
sociated with solving an inverse problem. For DLS (or 
XPCS), there are well known algorithms such as CON- 
TIN [531 for backing out P(a) via the distribution of 
decay times in the ISF. Or, less ambitiously, cumulant 
analysis [58j can be used to extract a. Form factors from 
static scattering are seldom inverted directly to yield 
P(a). Instead, one assumes, say, a Gaussian form, and 
the scattering profile from Mie theory is fitted to obtain 
(a) and a. The effect of small polydispersities (a few %) 
on the ISF and form factor can be treated analytically 
|59j . and becomes independent of the form of P(a) as 
a — > 0. The resulting expressions can be used to fit dy- 
namic or static scattering data to yield rather accurate 
values of a. 



V. MEASURING VOLUME FRACTION 



We now turn to describe and evaluate a number of 
methods for determining the volume fraction of model 
colloids. 
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A. Measuring mass and density 

The method used by Pusey and van Megen to deter- 
mine m their classic work on hard-sphere colloid phase 
behaviour Q and described subsequently in detail in a 
symposium paper [17f remains conceptually the simplest. 
They dried a suspension of known total (or 'wet') mass 
to determine the mass of dry particles, and converted the 
resulting mass fraction into using literature values of 
the densities of the solvents and of (dry) PMMA. There 
are multiple assumptions behind this procedure that lead 
to systematic uncertainties. In particular, this procedure 
assumes that the properties of dry and wet particles are 
the same, which, due to solvent absorption and solvation 
of the 'hairs', is unlikely to be true. Thus, many have 
subsequently proceeded differently. If has been sep- 
arately determined for one sample using other methods 
(e.g. at 4>f or 0m, see below), the the ratio of mass to 
volume fractions can be used to calibrate other samples. 
But the exact relationship between these two quantities 
is not a direct proportionality, and involves (unknown) 
ratios of the properties of wet and dry particles. 

A somewhat more involved procedure is in principle 
less problematic [60j . First, one determines the hydrody- 
namic radius an from dynamic light scattering. Then 
the sedimentation velocity, v s , of a dilute suspension is 
determined by analytic centrifugation to obtain the den- 
sity difference between the (wet) particles and the sol- 
vent, Ap = pj et - p s : v s = 2Apga 2 H /9r) , where 770 
is the solvent viscosity (separately measured) and g is 
the gravitational acceleration, although the assumption 
that the non-slip boundary condition holds at the 'hairy' 
particle surface may not be strictly valid (E. Sloutskin, 
personal communication). Since liquid densities can be 
determined very accurately using pyenometry or other 
densitometric methods if the temperature is controlled, 
we can measure p s and the density of an arbitrary sus- 
pension, p, from which its can be determined using 

p = p s + 0Ap. (5) 



B. Measuring phase behaviour 

A popular method of calibrating relies on the known 
phase behaviour of hard spheres. In particular, in the 
region 4>p = 0.494 < < m = 0.545, hard spheres 
show coexistence of fluid at 4>p and crystals at 0m- The 
fraction of crystals, %, increases linearly from to 100% 
over the interval. Measuring x f° r a sample within the 
coexistence region then gives its 0. To determine \ ac- 
curately, one needs to take into account the compression 
of the crystalline sediment by its own weight [6 lj. 

The main uncertainty associated with using phase be- 
haviour to calibrate is the effect of polydispersity. All 
calculations and simulations to date agree that finite a 
increases <j>p and 0m- Thus, e.g., in the 'moment free 
energy' calculations shown in Fig. @J 0f = 0.5074 and 



0m = 0.5540 at a = 5%, the latter being a representa- 
tive value of a typical preparation of PMMA colloids. To 
date there has been no independent experimental check 
on such theoretical predictions, one of the main issues 
being the measurement of in polydisperse colloids! 

Nevertheless, these results may throw light on one of 
the puzzles remai ning from the original work of Pusey 
and van Megen 0, Il7l | who found that if they assumed 
a freezing point of <pp = 0.494, their measured melting 
point was 0m = 0.535. The ratio 7 of these two values, 
which characterises the width of the coexistence gap, is 
7PvM = 1-083. For monodisperse colloids, 7 CT= o = 1-103, 
while calculations [HJ for a = 5% gives 7^=5% = 1.092. It 
is therefore possible that the narrowing of the coexistence 
gap observed by Pusey and van Megen is largely due to 
polydispersity Note, however, that the phase diagram 
likely depends on the whole P{a) and not just a. 

An additional source of uncertainty is residual 
charge [6^. [63| so that, hard-sphere phase behaviour no 
longer obtains: crystallization is expected at lower vol- 
ume fractions (for PMMA, see, e.g.[64|). In these cases 
the phase behaviour cannot be matched to that of hard 
spheres at all. However, by adding salt, the charges can 
be screened and hard sphere behaviour recovered to an 
extent (e.g., for PMMA, seej6|). 

We mention that confocal microscopy of a sample in 
the fluid-crystal coexistence region can be used to deduce 
a value for a e s by assuming particular values for (£>p and 
0m [HI, subject to all the above-mentioned caveats and 
uncertainties. 



C. Centrifugation and sedimentation 

Perhaps the quickest way to obtain samples with ap- 
proximately calibrated is by ccntrifuging to obtain a 
sediment that one assumes to be at 'random close pack- 
ing' (RCP), and therefore some known 0rcp, which can 
then be redispersed with fixed volumes of solvent to give 
samples at lower concentrations. The method can be 
applied even with charged particles, since hard centrifu- 
gation can reduce even such particles to a mutually- 
touching amorphous state [43|. 

The main problem with this method is that the theo- 
retical status of RCP is still debated, with different simu- 
lation algorithms giving different results [65|, [66j|. Exper- 
imentally, different regions of the centrifuged sediment 
have somewhat different concentrations (0 = 0.60 — 0.64 
in silica colloids [67j), and little is known about the 
almost-certain dependence of sediment structure on cen- 
trifugation protocol. Moreover, the spun-down sediment 
is inevitably compressed, and will expand with time af- 
ter the cessation of centrifugation, which introduces an 
extra degree of uncertainty. Finally, the dependence on 
polydispersity is poorly known [66l l68l. |69| . 

But centrifugation is convenient, and if the protocol is 
kept constant, it can be used to produce a series of sam- 
ples with highly accurate normalized concentrations, viz., 
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4>/4>sed, where se d is the volume fraction of the sediment. 

Under this heading, we may mention that particles 
with small enough gravitational Peclet number [43l |70| 
(cither by virtue of near density matching or by virtue 
of being small) and low enough polydispersity will sedi- 
ment slowly under gravity to form sedimentary crystals 
consisting of more or less randomly-stacked hexagonal 
close packed (rhcp) layers of particles. If the particles are 
monodisperse hard spheres, then </> r h cp = 7r/vl8 ps 0.74 
in this sediment. Again, however, the (largely unknown) 
effect of polydispersity as well as any changes due to 
charges need to be taken into account. 



D. Confocal microscopy and particle counting 

Confocal microscopy can be used to locate the posi- 
tion of thousands of particles in a suspension. Thus, if 
the particle radius, a, is known, then counting TV par- 
ticle in an imaging volume V will yield <j) directly using 
Eqn. [1] Occasional particle mis-identification or missing 
a particle all together by the software give rise to erro- 
neous <j), so that it is important to cross-check particle 
positions identified against raw images. In particular, 
particles near the edge of images are often mis-identified, 
so that in practice a sub-volume only is considered. Fi- 
nally, uncertainties in a are magnified 3-fold or more in 
calculating <f>. This latter uncertainty is compounded by 
the issue of which of the possible radii (Section llll Al) one 
should use. 



E. X-ray transmission 

The intensity of X rays transmitted by a sample is 
given by It = Ia& ^ x , where Iq is the incident inten- 
sity, [i and x are the attenuation coefficient and thick- 
ness of the sample. In the case of a colloidal suspension, 
fi = (1 — 4>)fjis + <^cA*pi where <pc is the volume fraction of 
particle cores, and fi s and \i v are the attenuation coeffi- 
cients of the solvent and particles. The negligible amount 
of electron density represented by stcrically-stabilizing 
'hairs' means that they hardly contribute to the beam 
attenuation. X ray transmission can therefore be used to 
determine <p directly for model colloids such as charge- 
stabilised polystyrene [7l[ or silica[72j, but only the core 
volume fraction for sterically-stabiliscd particles. 



F. Measuring ^-dependent properties 

The ^-dependence of a number of material properties 
of hard-sphere suspensions are known either from ana- 
lytic theory or highly-accurate simulations. In principle, 
therefore, measuring these properties can be used to de- 
termine <f>. Here we review three: viscosity, diffusivity 
and structure factor. 



Einstein predicted that in the limit </>—»• 0, the vis- 
cosity of a hard-sphere suspension is given by /?7n = 
1 + (5/2)</j, with rjo being the viscosity of the solvent [34| . 
Thus, in principle, measuring r)((f>) is a method for deter- 
mining <f> (e.g.[73[). While suspensions in general shear 
thin, this should not be a problem in the very dilute 
limit. But temperature control is important, since 770 is 
temperature sensitive (cf. Section fill Bp . 

The problems associated with this method have been 
detailed beforejll]. In essence, very low <f>, certainly 
< 0.02, must be reached for the Einstein result to be 
valid; otherwise, second[75|, [ZfjJ and higher order term in 
this 'virial' expansion needs to be taken into account. In 
the case cited[73|], using the Einstein relation at <fi w 3% 
leads to an error in 4> of ~ 7%[74|. The difficulty, of 
course, is that in the limit 4> — >• 0, very accurate viscom- 
etry is needed to distinguish the dilute suspension from 
pure solvent. Using the Einstein relation to calibrate 4> 
in suspensions that are too concentrated for the relation 
to be valid accounts for some of the spread in literature 
values of ry(^), the viscosity of the most concentrated 
stable fluid state of hard spheres. Interestingly, determin- 
ing <f> using the Einstein relation is strictly independent 
of polydispersity: in the dilute limit, each particle con- 
tributes by an additive amount that is proportional to its 
volume. 

Instead of measuring (</>), one could determine the 
single-particle diffusion coefficient as a function of cf>. 
Thus, El Masri et aZ. measured the short-time self dif- 
fusion coefficient as a function of volume fraction, ((f) . 
The difficulty is that there are at least two different pre- 
dictions for this behaviour [z||, [7^] which leads to a 7% 
absolute uncertainty in </>. 

Lastly, we have already mentioned ( Section llll Bp that 
analytical expressions for the static structure factor, 
S(q), of hard spheres are available. In particular, the 
closed-form expression from the Percus-Yevick (PY) ap- 
proximation [35J fits simulation data closely, provided 
that the empirical Vcrlet-Weis correction to the volume 
fraction[8(| is applied, i.e. the PY structure factor for 
volume fraction 0' is used for an experimental sample at 
<j>: <t>' = <fi— 4> 2 /16. Thus, fitting measured S(q) can yield 
a measure of <p, provided that the particles can be treated 
as hard spheres. Again, caution about residual charges 
applies. Alternatively, g(r) determined from confocal mi- 
croscopy can be fitted to the PY form or to simulation 
data[81| to give <\>. 



G. Deceptive samples 

Finally, we explain how using an accurately calibrated 
'stock colloid' may still lead to errors in the 4> of samples. 

First, we have already mentioned a number of times the 
issue of swelling. If particles used for calibrating volume 
fraction are still in the process of swelling due to solvent 
absorption, then samples prepared subsequently will have 
a higher </> than the earlier calibration would suggest. 
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Secondly, preparing samples almost invariably involves 
transferring suspension from one container (e.g. a bot- 
tle of stock) to another (e.g. a capillary for microscopy) 
using (typically) a pipette or a syringe. Apart from dif- 
ficulties caused by very high viscosities [22[ and shear 
thickening [sH, [83|], there is the problem of jamming of 
the particles as the suspension enters a constriction [84j . 
which leads to a 'self filtration' effect. Particles jammed 
at (say) the entrance to a pipette prevents other particles 
from entering, but solvent continues to flow, so that the 
sample inside the pipette has a lower <f> than the bulk 
suspension that we hope to transfer. Thus, a sample 
loaded for confocal microscopy may be more dilute than 
one expects. 

Thirdly, except for very well density-matched sam- 
ples at a temperature accurately remaining at the tem- 
perature at which the density matching was originally 
achieved, suspension inevitably sediment (or cream) with 
time at all except </>r,cp or r hcp- This will lead to con- 
centration gradients. Indeed, such gradients can be de- 
liberated exploited (43l lull l85l - [87| . e.g. to determine 
equations of state. But in other cases, concentration gra- 
dients lead to unintended local deviations from the av- 
erage 4> at which the sample as a whole was originally 
prepared. 

Since one of the most important uses of hard-sphere 
colloids is as a model to study dynamical arrest fl5l [88L |89| 
and associated properties such as aging [90ll9lj |. any of the 
above three sources of unintended changes in <p will have 
severe consequences: all suspension properties change 
very rapidly with <j> at and above the glass transition 
(</> > 0.58). 



H. Summary: relative vs. absolute (j> 



The most important message from the preceding crit- 
ical review is that the statistical errors involved in de- 
termining (j) in the competent use of any of the above 
methods can almost certainly be brought below the sys- 
tematic errors involved. Thus, for example, if one uses 
confocal microscopy to count particles and thus deter- 
mine <f>, the major source of uncertainty is likely to be the 
input radius. Thus, it is perfectly possible to produce a 
series of samples with relative uncertainty in of 1 part 
in 10 . However, our collective experience in using many 
of these methods suggests that the systematic uncertain- 
ties are unlikely to be below 3-6%. Far from being a 
small error, such uncertainties can have dramatic effects. 
Thus, e.g., the viscosity of a hard-sphere suspension [74j 
grows by a factor of 2 when <j> increases from 0.47 to 0.49; 
and the simulated crystal nucleation rate [94| near <frp can 
change by 10 orders of magnitude for a 1% change in the 
absolute value of 6. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 



FIG. 5: (a) Raw data: (■): The normalized long-time self 
diffusion coefficient of PMMA colloids as a function of volume 
fraction, Do/D f (</>)[§2], with the volume fraction multiplied 
by 1.022. (A): The normalized low-shear viscosity of PMMA 
colloids, r)(4>)/r]o, as a function of volume fr action [22|. (♦): 
the volume fraction of the viscosity data set being multiplied 
by 0.965. (b) Violation of Stokes-Einstein-Sutherland rela- 
tion. The factor f3, Eqn. [7J from the raw data (♦), from 
theory[93| (B) and after the <j> values of the viscosity data set 
having been multiplied by 0.965 (a). 



VI. A CAUTIONARY TALE 

In this section we give a case study to illustrate how 
important it is to be critical about experimental <f> values 
by analyzing two published data sets. These data sets 
give as a function of <fi the long-time self diffusion coef- 
ficient, D^(ip) [92|, and the low-shear viscosity, [2^]. 
of PMMA colloids. At the time of publication, the dif- 
fusion data was the best and most complete available, 
and the viscosity data remain one of the most complete 
to date. The two groups came to quite different conclu- 
sions about dynamical divergence at high <f>. van Megen 
and his collaborators concluded that Dg(cf)) diverged at 
(f> = (f) g w 0.58, in a manner consistent with that pre- 
dicted by mode coupling theory for an ideal glass transi- 
tion. Chaikin and his collaborators, however, concluded 
from their rj((f>) that there was no glass transition at 
« 0.58; instead, they suggested that rj((j)) diverged at 
RCP, 4> = 0.64, according to a Volgel-Fulcher law. This 
controversy is ongoing (see, e.g., [95|,[96|]). We do not to 
enter into this discussion here; instead, we use older data 
to illustrate many of the issues concerned with measuring 
<f> and using experimental data sets. These issues are, of 
course, pertinent for the ongoing discussion. 

The first thing to notice about these two experiments 
is that the reported volume fractions cannot be com- 
pared directly. Both sets of authors relied on measuring 
phase behaviour to calibrate <fi (cf. Section lVB[) . but one 
set of authors took into account polydispcrsity, and one 
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did not. van Megen and co-workers used 4>p — 0.494, 
the value for monodisperse hard spheres, to determine 
0, but cautioned that their particles had a polydispcr- 
sity of a = 5%. The colloids used by Chaikin and his 
co-workers also had the same <r, and they used the sim- 
ulation data of Bolhuis and Kofke [53[ to move freezing 
to <j>p = 0.505 for this polydispersity. Interestingly, the 
latest analytic calculations agree closely: Wilding and 
SollichfH give <\> F = 0.5074 at a = 5%, Fig. ffl Thus, 
we multiply the (f> value of the van Megen data set by a 
factor of 0.505/0.494 = 1.022 to make it consistent with 
the Chaikin <f> values. The resulting data are shown in 
Fig. El The measurements have been normalised, D^((f) 
by the single-particle diffusivity, D , and r]{<f>) by the sol- 
vent viscosity, rjo . The normalized viscosity diverges at 
higher (f> than the normalized (inverse) diffusivity. 

At (/)—)■ 0, the solvent viscosity and the single-particle 
diffusivity are related by the Stokcs-Einstcin-Suthcrland 
relation (SESR): Dq = k B T '/Qirrjoa. At finite cj>, there is 
no a priori reason that a generalized SESR should hold 
for any of the many diffusion coefficients that can be 
defined. So we write 



£ S L (</>) = 



k B T 
67i7?(0)a 



x (3 



(6) 



where j3 is a numerical factor that can be restated as 
= Dfa) x 77(0) 



Dn 



(7) 



We plot in Fig. (5Jb) (diamonds) the /3((f>) implied by the 
data sets in Fig. Eta). In so far as ft ^ 1, the SESR is 
violated. 

Violation of the SESR is widely known for glass- 
forming systems near the glass transition. In all experi- 
mental cases known (see e.g. [13), /3 > 1, i.e. the particles 
diffuse somewhat faster than the viscosity allows accord- 
ing to the SESR. The fact that /3 drops very substantially 
below unity at <f> > 0.4 in Fig.[SJb) is therefore surprising, 
and merits further analysis. 

To proceed, we turn to the work of Banchio et al. [93], 
who have calculated various diffusivities and viscosities of 
hard sphere suspensions within a mode-coupling frame- 
work, and have shown that their results compared well 
with multiple experimental data sets. Their calculations 
predict that (3 as defined in Eqn. [7] hovers just below 
unity in the range < <fi < 0.50. Fig. [5jb) shows that 
the experimental (3(4>) from the data plotted in Fig. [5Ja) 
(diamonds) essentially agrees with theory (squares) up 
to 4> = 0.35, but start to diverge thereafter. 

Since we conclude that absolute values of </> are unlikely 
to be accurate to better than 3 — 6%, it is interesting to 
note that multiplying the volume fractions in the viscos- 
ity data set by a factor of 0.965 overlaps the two nor- 
malized data sets, Fig. [SJa). Not surprisingly, then, this 
rcnormalization of <j) also brings very substantially bet- 
ter agreement in /3(0) in the whole range of </> covered by 



theory [93 (triangles, Fig. Efb)). Assuming that there is 
a glass transition at (f> g w 0.58, then this rcnormalization 
of (j> also brings the direction of SESR violation in the 
vicinity of 4> g in line with all other known glass formers, 
viz., (3 > 1. Thus, the supposed disagreement between 
the two data sets is well within the range of expected 
uncertainties in the absolute determination of d>. 



VII. CONCLUSION 

Hard sphere colloids are now part of the accepted 'tool 
kit' of experimental statistical mechanics. What we aim 
to do in this critical review is to counsel caution in com- 
paring data from experiments against theory or simula- 
tions, because there are substantial, and probably irre- 
ducible, systematic errors in determining suspension vol- 
ume fraction. This situation calls for at least three re- 
sponses. First, experimentalists need to take the cue from 
the pioneering work of Pusey and van MegenfT?) and al- 
ways report exactly how they arrive at their quoted <f> 
values, and discuss likely sources particularly of system- 
atic errors. Secondly, experimental data sets need to be 
compared vigilantly against each other to reveal possible 
discrepancies. Finally, theorists and simulators seeking 
experimental confirmation of their results should not be 
too easily satisfied with apparent agreement, at least not 
until in-depth inquiry into the systematic uncertainties 
in cj) has been carried out. 

Finally, we note that while perfect hard spheres are 
indeed characterized by a single thermodynamic variable 
<t>, real particles are never truly hard. Some softness in 
sterically-stabilized particles necessarily comes from com- 
pressible 'hairs', but this becomes less significant as a 
increase. However, it is becoming clear that for larger 
(a > 0.5/i/m) PMMA particles, a certain degree of charg- 
ing is inevitable [62], |63|, |8l|, which cannot be entirely 
screened by salt (due to limited solubility in organic sol- 
vents). Such softness means that accurate measurement 
of <p alone is insufficient, and introduces further uncer- 
tainties. 
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